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HIGHLIGHTS 


• We study the real-time analysis of chemical compositions of corn stover using NIRS. 

• NIRS can be used for real-time analysis of these chemical compositions. 

• These compositions could be predicted with varying accuracy using NIRS. 

• This method has potential for industrial applications. 
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On-line analysis of proximates (moisture, ash, volatile matter, and fixed carbon) and lignocellulose 
components (cellulose, semi-cellulose, and lignin) of coarse-crushed corn stover was investigated using 
near-infrared spectroscopy (NIRS). Compared with traditional technique, on-line measurement is com¬ 
plex and challenging but provides real-time analysis. The spectrometer operated over a conveyor belt 
that moved at 20 cm/s, and the distance between the spectrometer scanning window and the surface 
of the sample was 100 mm. Corn stover samples (n = 217) were collected from three provinces in China 
and used to develop the models. Samples were crushed to <50 mm before analysis. After optimized pre¬ 
treatment, all the NIRS models were developed using partial least squares (PLS) method. The relative 
standard deviations (RSD) of the models for moisture, ash, volatile matter, fixed carbon, cellulose, 
semi-cellulose, and lignin were 9.04%, 9.99%, 1.31%, 6.00%, 3.87%, 6.80%, and 5.06%, respectively. On-line 
analysis of proximates and lignocellulose components of corn stover was possible using this NIRS system. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stover is the stalk residue left after crop production, and can be 
used as a renewable and environmental friendly alternative source 
of energy [1]. Corn stover is a major field crop residue in China, 
accounting for 46% of all stover and with annual production of over 
172 million tons [2 . It can be used to produce bioethanol, pyro- 
lyzed for fuel gas, liquefied for biodiesel, torrefied for upgraded 
solid fuel, converted into hydrochar under hydrothermal carbon¬ 
ization (HTC) process, and combusted to produce electricity and 
thermal energy [3-7]. However, its suitability for energy use can 
be affected by the crop location, plant species, collection method, 
and storage time, which can all cause variation in the corn stover 
composition, including proximates (moisture, ash, volatile matter, 
and fixed carbon contents) and lignocellulose components (cellu¬ 
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lose, semi-cellulose, and lignin contents). Therefore, to optimize 
energy production, a method for determining the corn stover 
composition before use is required. For example, the moisture, cel¬ 
lulose, semi-cellulose, and lignin contents of stover need to be 
known for pretreatment and fermentation for production of 
bioethanol and biodiesel. In addition, the moisture content, ash 
content, volatile matter, and fixed carbon are required for reactor 
design and process technology development for pyrolysis and com¬ 
bustion of biomass. 

Unfortunately, traditional analysis of proximates and lignocel¬ 
lulose components is time consuming, labor intensive, and expen¬ 
sive, and on-line analysis is difficult. By comparison, near-infrared 
spectroscopy (NIRS) is fast, efficient, environmentally friendly, and 
nondestructive. Near-infrared spectra reflect the overtones and 
combination of fundamental vibrations of functional groups such 
as NH, OH, and CH [8]. As a rapid alternative to traditional tech¬ 
niques, NIRS has been successfully used for analysis of proximates 
and lignocellulose components analysis of biomass [9-15]. 
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However, most of these earlier studies have used laboratory 
spectrometers. For industrial applications, there is an urgent need 
for on-line detection in real-time, which provides the ability to 
detect the raw material, process changes, and end product at the 
same time. On-line detection also improves the analysis efficiency 
and provides prompt feedback for process control. These benefits 
can decrease the manpower requirements, and increase the safety 
and integrity of the analysis. However, collection of spectra is dif¬ 
ficult with a moving sample [16]. In addition, noise can increase 
because of the spectral collection method and mechanical vibra¬ 
tions of the collection system. Consequently, on-line detection of 
proximates and lignocellulose components of corn stover using 
NIRS is challenging. What is more, for most industry application, 
like pyrolysis for fuel gas and producing bioethanol, the corn stover 
was normally use harmer mill and use screen to troll the particle 
size. However, utilizing of biomass for combusting to produce 
electricity and thermal energy, the corn stover are simply crushed 
(size range being 10 mm, 20 mm or even larger) [17]. And the corn 
stover was often coarsely crushed before being transported to the 
plant with different kinds of mill. 

Based on the above consideration, the potential of NIRS for on¬ 
line determination of proximates and lignocellulose components of 
coarse-crushed corn stover was investigated in this study. 

2. Material and methods 

2.2. Samples 

Two hundred and seventeen corn stover samples of different 
varieties were collected from three provinces (seven locations 
per province) in China (Hebei, Shanxi, and Shandong) in 2012. After 
air drying, the samples were crushed in a tooth mill (without a 
sieve) until most of the fragments were shorter than 50 mm 
(Fig. 1A). Furthermore, according to the results of the experiment 
(sieving method), 98.17% (W/W) of the sample can pass through 


a sieve having a mesh aperture of 20 mm (Fig. IB) and 69.04% 
(W/W) of the sample can pass through a sieve having a mesh 
aperture of 10 mm (Fig. 1C). The crushed samples were used for 
spectral analysis. For chemical analysis, some of the samples were 
ground to pass through a 1 mm screen (Fig. ID). Prepared samples 
were stored in a dry and cool room before analysis. 

2.2. Reference analysis 

To establish reference values for the parameters, the samples 
were first analyzed according to standard methods (Table 1 ). The 
fixed carbon was calculated by subtracting the sum of the percent¬ 
age moisture, ash, and volatile matter from 100. All the values were 
the average of two measurements and are reported in g/kg. All the 
values were reported on as-received basis. 

2.3. On-line NIRS spectra collection 

The on-line measurement system consisted of a NIR spectrom¬ 
eter, a conveyor belt, a sample holder, and a computer (Fig. 2). 
Matrix-F (Bruker Daltonik GmbH, Bremen German) was used to 
collect the spectra. A fiber-coupled sensor head was placed on 
the top of the sample holder, which sat over the conveyor belt, 
and the entire machine was placed on a table with the computer. 
To collect sample spectra, the sensor was passed through a hole 
(ol45 mm, 100 mm above the conveyor belt) cut into the sample 
holder. Spectra were collected in diffuse reflection mode between 
12,000 and 4000 cm -1 at a resolution of 8 cm -1 . An internal refer¬ 
ence (golden board) was used. We had studied the acquisition 
parameters previously [24]. In the present study, the parameters 
were determined using the results. To ensure the distance between 
the sample and the scanning window was consistent, each sample 
was placed on a steel box (700 (length) x 160 (width) x 50 
(height) mm). The steel box was moved forward by the conveyor 
belt. The distance between the scanning window and the surface 



Fig. 1 . Corn stover samples prepared for spectral and chemical analysis. (A) Coarse-crushed corn stover, (B) undersize (20 mm) and oversize (10 mm), (C) undersize (10 mm) 
and (D) ground corn stover. 
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Table 1 

Chemical analysis methods. 


Parameters Method or references 


Instrument 


Preparation of 
sample 
Moisture 
Ash 

Volatile matter 

Cellulose 

Semi-cellulose 

Lignin 


ASTM E1757-01 standard practice for preparation of biomass for compositional analysis [18] Grinding machine 


ASTM El 756-08 standard test method for determination of total solids in biomass [19] 

ASTM El 755-01 (2007) standard test method for ash in biomass [20] 

ASTM E872-82 (2006) standard test method for volatile matter in the analysis of particulate 
wood fuels [21] 

NREL/TP-510-42618 determination of structural carbohydrates and lignin in biomass [22] 

ASTM El 721-01 standard test method for determination of acid-insoluble residue in biomass 

[23] 


Drying oven 
Muffle furnace 
Muffle furnace 

High performance liquid chromatography 
(Agilent, USA) 

Muffle furnace 



Fig. 2. Spectra collection system. 


of the samples was 100 mm and the speed of the belt was 20 cm/s. 
When the sample box moved underneath the scanning window, 
scanning began, and it ended when the required number of scans 
were acquired. Each spectrum is the average of 30 scans. If the 
required number of scans had not been acquired by the time the 
sample box moved out from under the scanning window, the pro¬ 
cess was repeated as many times as required. In these cases, the 
average spectrum was used as the final spectrum. Each sample 
was reloaded and measured three times, and the three spectra 
were averaged. 

2.4. Development and analysis ofNIR models 

All NIR models were developed using the partial least squares 
(PLS) Toolbox of MATLAB R2008a (Mathworks, Natick, MA, USA). 
At first, samples were divided into calibration and validation sets 
by taking every fourth sample from the sample set ascending in 
reference data into validation set, and the remaining samples 
comprised the calibration set. Then the spectra were pretreated 
using an autoscale, first derivative or second derivative, smoothing 


(SavGol), standard normalized variate (SNV), multiplicative scatter 
correction (MSC), or a combination of these manipulations. For 
each original and pretreated spectra, a PLS regression was 
calculated to determine the quantitative relationship between 
the spectral variables and the stover composition. 

The final model was chose according to Rc 2 , the root mean 
square error of calibration (RMSEC), Rp 2 , and the root mean square 
error of prediction (RMSEP). The accuracy of the model was evalu¬ 
ated using R 2 , RMSEC, RMSEP, and the relative standard deviation 
(RSD). The RSD was calculated according to formula (1): 

RSD = RMSEP/Mean x 100%, (1) 

where Mean is the average of the chemical analysis value in the val¬ 
idation set. 

3. Results and discussion 

3.1. Reference parameters 

The range, mean, and standard deviation for the reference 
analysis are presented in Table 2. Several samples were excluded 
from the sample set because of redundancy or they were outliers. 
There were 5, 3, 7, 5, 11, and 11 samples excluded for moisture, 
volatile matter, fixed carbon, cellulose, semi-cellulose, and lignin, 
respectively. Histograms of the parameters (Fig. 3) showed all 
the parameters were normally distributed. For the proximates, 
volatile matter was the main component and accounted for 
approximately 75% of the total weight of the sample. This is consis¬ 
tent with earlier results [25]. Generally, corn stover contains more 
volatile matter and less fixed carbon than coal. Consequently, it has 
a lower ignition temperature and lower heating value than coal. 
Cellulose, semi-cellulose, and lignin combined accounted for 
approximately 66% of the total weight of the sample, which agreed 
with reported values [10,12]. 

3.2. Spectra 

Spectra of the samples are presented in Fig. 4A. Because of noise 
in the 4000-4440 cm -1 range and the lack of useful information in 


Table 2 

Statistical analysis of the calibration and prediction set. 


Parameters 

Calibration 



Validation 



Range 

Mean 

S.D. a 

Range 

Mean 

S.D. 

Moisture (g/kg) 

24.85-70.80 

41.11 

7.23 

25.97-70.02 

41.64 

7.41 

Ash (g/kg) 

42.06-144.55 

75.25 

18.81 

42.99-125.27 

75.56 

18.70 

Volatile matter (g/kg) 

674.76-797.32 

754.15 

19.23 

714.91-798.59 

754.15 

19.23 

Fixed carbon (g/kg) 

81.47-167.74 

130.26 

13.90 

100.54-158.63 

130.59 

13.44 

Cellulose (g/kg) 

320.06-474.80 

394.54 

31.02 

327.75-465.01 

395.20 

30.51 

Semi-cellulose (g/kg) 

103.46-175.41 

139.90 

14.99 

103.96-170.75 

139.29 

15.20 

Lignin (g/kg) 

188.38-295.83 

233.30 

20.39 

194.14-280.12 

231.78 

19.01 


a S.D. represents standard deviation. 
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Moisture (g/kg) Ash (g/kg) 



670683696709722735748761774787800 100107114121128135142149156163170 

Volatile matter (g/kg) Fixed carbon (g/kg) 




Cellulose (g/kg) Semi-cellulose 



190201212223234245256267278289300 

Lignin (g/kg) 

Fig. 3. Frequency histograms of reference parameters: (A), (B), (C), (D), (E), (F) and (G). 

the 8000-12,000 cm -1 range (not shown), the spectra only show of water (6897 cm -1 ), O—H stretch first overtone of cellulose 
4440-8000 cm -1 . There were three main peaks between 4440 (6711 cm -1 ), C—H combination band of lignin (6897 cm -1 ), and 

and 8000 cm -1 , and these arose from cellulose, semi-cellulose, lig- N—H stretch first overtone of protein (6623 cm -1 ). The peak at 

nin, and moisture. The wide peak located between 6500 and around 5200 cm -1 was caused by the O—H stretch/HOH deforma- 

7000 cm -1 could be attributed to the O—H stretch first overtone tion combination of cellulose. Finally, the peak at around 
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(A) 0.9 



8000 7000 0000 fO>D 4000 

Wavenumbetfcm' 1 


Samples/Scores Plot 



-50 0 50 100 


Scores on PC 1 (91.65%) 


Fig. 4. Spectra and PCA analysis of the spectra of samples. (A) Raw spectra of the 
samples and (B) principal components (PCI and PC2) of the spectra. 

4800 cm -1 was linked to the asymmetrical C—O—O stretch third 
overtone (6897 cm -1 ) of cellulose, and the symmetrical N—H 
stretch/amide combination (4866 cm -1 ) and N—H bend second 
overtone or N—H bend/N—H stretch of protein (4854 cm -1 ). 

Principal component analysis (PCA) can be used to divide sam¬ 
ples into several groups with similar characteristics and to identify 
possible outliers [26]. PCA was applied after autoscale pretreat¬ 
ment was used to standardize the spectral data. According to the 
PCA plot (Fig. 4B), the first two PCs explained 98.39% of the total 
variation in the spectral data. The sample distribution was very 
wide, which meant the samples were quite representative. 


3.3. On-line NIRS models 

After appropriate pretreatment (Table 3) of the spectra, the best 
prediction models for the proximates and lignocellulose 

Table 3 

Results for the calibration and prediction models. 


components of corn stover were established using PLS_Toolbox. 
The statistical results of the models are shown in Table 3, and scat¬ 
ter plots of predicted versus measured values are shown in Fig. 5. 
The variable importance in the projection scores of wavenumbers 
for the models are presented in Fig. 6. These values indicate the 
wavenumber regions in the spectra that are most sensitive to 
changes in the stover composition. This information is important 
for establishing robust models and developing spectrometers for 
specific applications. 

3.3.1. Prediction of proximates (moisture, ash, volatile matter, and 
fixed carbon contents) 

For moisture, the Rc 2 , RMSEC, RMSEP, and RSD were 0.81, 
3.18 g/kg, 3.76 g/kg, and 9.04%, respectively. It was possible to 
use the model to predict the moisture of the sample. Important 
spectral regions for the prediction were identified (Fig. 6A). Absor¬ 
bance bands in the 5000-5600 cm -1 region corresponded to the 
O—H bend second overtone (5155 cm -1 ) and O—H combination 
band (5587 cm -1 ) of water. The bands in the 6600-7200 cm -1 
region corresponded to the O—H stretch first overtone 
(6711 cm -1 ) and O—H stretch first overtone (6897 cm -1 ) (Fig. 6A). 

The model could also be used to predict the ash content 
(Rc 2 = 0.85, RMSEC = 7.21 g/kg, RMSEP = 7.55 g/kg and 
RSD = 9.99%). Simple inorganic compounds will not absorb in the 
NIR region, but inorganic complexes will [27]. Compared with 
organic substances, the spectral bands for inorganic substances 
are broader and fewer in number. The carbonate ion typically 
absorbs in the NIR region, and in the present study bands were 
observed at around 4630 cm -1 (Fig. 6B). 

For volatile matter, three peaks were observed that corre¬ 
sponded to the O—H first overtone (7092 cm -1 ), C—H stretch first 
overtone (5797 cm -1 ), and O—H bend/C—O stretch combination/ 
asymmetrical C—O—O stretch third overtone (4762 cm -1 ) 
(Fig. 6C). The model could be used to predict the volatile matter 
content, with a RMSEP of 9.92 g/kg and RSD of 1.31%. 

For prediction of the fixed carbon content, many peaks were 
useful. These were associated with the C—H stretch/C=0 stretch 
combination (4545 cm -1 ), C=0 stretch second overtone 

(4926 cm -1 ), C—H stretch first overtone (5935 cm -1 ), and C=0 
stretch third overtone (6897 cm -1 ) (Fig. 6D). The model could be 
used to predict the fixed carbon content (Rc 2 = 0.76, 
RMSEC = 6.87 g/kg, RMSEP = 7.83 g/kg, and RSD = 6.00%). The accu¬ 
racy of the model for this prediction was not good, perhaps 
because the fixed carbon was calculated from three other values 
and would have an accumulated error. 

3.3.2. Prediction of lignocellulose components (cellulose, semi¬ 
cellulose, and lignin contents) 

For cellulose, the model could be used for prediction with Rc 2 , 
RMSEP, and RSD of 0.82, 15.28 g/kg, and 3.87%, respectively. The 
wavelengths for the prediction of cellulose corresponded to C—H 
stretch first overtone (5797 cm -1 ), C—H stretch first overtone 


Parameters Calibration Validation 



Pretreatment 

PLS factors 

Rc 2 

RMSEC (g/kg) 

Rp 2 

RMSEP (g/kg) 

RSD (%) 

Moisture 

Smoothing (0, 3) a , baseline (weighted least squares), autoscale 

6 

0.81 

3.18 

0.75 

3.76 

9.04 

Ash 

1st Derivative (2, 21 ) b , autoscale 

11 

0.85 

7.21 

0.84 

7.55 

9.99 

Volatile matter 

SNV, smoothing (0, 9), 1st derivative (2, 25) 

10 

0.67 

10.94 

0.68 

9.92 

1.31 

Fixed carbon 

Smoothing (3, 3), 1st derivative (2, 3), autoscale 

7 

0.76 

6.78 

0.66 

7.83 

6.00 

Cellulose 

Smoothing (0, 3), 1st derivative (2, 3), autoscale 

9 

0.82 

13.12 

0.77 

15.28 

3.87 

Semi-cellulose 

SNV, smoothing (0, 3), 1st derivative (2, 7), autoscale 

8 

0.66 

8.72 

0.62 

9.47 

6.80 

Lignin 

Smoothing (3, 5), 1st derivative (2, 5), autoscale 

7 

0.63 

12.31 

0.61 

11.73 

5.06 


a Smoothing (0, 3): “0” and “3” represent polynomial order and filter width, respectively. 
b 1st Derivative (2, 21): “2” and “21” represent polynomial order and filter width, respectively. 
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Measured (g/kg) Measured (g/kg) 



Measured (g/kg) Measured (g/kg) 



Measured (g/kg) 


Measured (g/kg) 



Measured (g/kg) 


Fig. 5. Scatter plot of predicted versus measured values: (A), (B), (C), (D), (E), (F) and (G). 
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Fig. 6. Variable importance in the projection scores of the wavenumbers for the models: (A), (B), (C), (D), (E), (F) and (G). 
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(5935 cm -1 ), C—H stretch first overtone (6173 cm -1 ), 0—H stretch 
first overtone (6494 cm -1 ), and O—H stretch first overtone 
(6711 cm -1 ) (Fig. 6E). 

The model for semi-cellulose could also be used for successful 
prediction, with Rc 2 , RMSEP, and RSD of 0.66, 9.47 g/kg, and 
6.80%, respectively. Only two peaks were identified for the predic¬ 
tion of semi-cellulose (Fig. 6F). These included a C—H stretch/C=0 
stretch combination of —CHO (4545 cm -1 ) and a C—H stretch first 
overtone/C—H stretch/HOH deformation combination (5618 cm -1 ) 
(Fig. 6F). 

The model for prediction of lignin content was also successful. 
The peaks identified for this prediction were a C—H stretch first 
overtone (5666 cm -1 ), C—H stretch first overtone of CH 2 
(5797 cm -1 ), C—H stretch first overtone of - CH 3 (5865 cm -1 ), and 
0—H stretch first overtone (6494 cm -1 ) (Fig. 6G). 

3.3.3. Comparison and analysis 

In general, the models in this study were not as accurate as 
those in earlier studies [5-11]. This could be attributed to the 
on-line detection, which increases the sources of error. First, the 
open sample box may bring external (ambient) light, vapors 
(steam, water vapor), and dust into the system [28], which will 
all decrease the accuracy of the model. Second, mechanical vibra¬ 
tions of the collection system will increase noise, shift spectral 
peaks, and change the optical density. Finally, the sample particle 
size (<50 mm) in this study may have been more variable than in 
the earlier studies. However, this particle size was chosen as it 
reflects the particle sizes that would be encountered in industrial 
use [17]. The particle size influences the optical path, which affects 
the absorbance and scattering coefficient of the sample to light. 
Overall, these sources of error would combine to reduce the predic¬ 
tion accuracy of the model. With the improvement of equipment 
and technological process, using sample with controlled size, the 
accuracy of the models will very likely be improved. 

4. Conclusions 

Successful NIRS on-line models for predicting proximates and 
lignocellulose components of coarse-crushed corn stover were 
established. Samples were ground to <50 mm and placed on a con¬ 
veyor belt moving at 20 cm/s, which moved them under a scanning 
window (distance to sample surface = 100 mm). The RSDs of the 
models for moisture, ash, volatile matter, fixed carbon, cellulose, 
semi-cellulose, and lignin contents were 9.04%, 9.99%, 1.31%, 
6.00%, 3.87%, 6.80%, and 5.06%, respectively. In future studies, the 
accuracy and robustness of the models should be improved before 
industrial application. 
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